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1 Introduction 



Several experimental collaborations at ep and pp colliders presented data on the differ- 
ential cross section d 2 a/dydp T for the inclusive production of D°, D + , and D+ mesons, 
A+ baryons, and their charge-conjugate counterparts. At DESY HERA, such data were 
collected by the ZEUS Collaboration [1,2] in low-Q 2 ep collisions, equivalent to photopro- 
duction, and by the HI Collaboration [3] in deep-inelastic ep scattering. At the Fermilab 
Tevatron, such data were taken by the CDFII Collaboration [4] in pp collisions. 

On the theoretical side, fragmentation functions (FF's) for the transitions c,b — > X c , 
where X c denotes a generic charmed hadron, are needed as nonperturbative inputs for the 
calculation of all the cross sections mentioned above. Such FF's are preferably constructed 
by using precise information from e + e~ — > X c + X via e + e~ annihilation at the Z-boson 
resonance, where X denotes the hadronic rest. In this process, two mechanisms contribute 
with similar rates: (i) Z — > cc decay followed by c — > X c (or c — > X c ) fragmentation; and 
(ii) Z — > bb decay followed by b — > X b (or b — > X b ) fragmentation and weak X b — > X c + X 
decay of the bottom-flavored hadron X b . The latter two-step process is usually treated 
as a one-step fragmentation process b — > X c . 

Using ALEPH [5] and OPAL [6] data on inclusive D* + production at the Z-boson 
resonance, we determined separate FF's for c — > D* + and b — > D* + in collaboration with 
Binnewies [7]. It is the purpose of this work to extract nonperturbative FF's for c,b — > 
D°,D + ,Df,A+ from the respective data samples collected by the OPAL Collaboration 
at LEP1 [8] using the same theoretical framework as in Ref. [7]. 

The work in Ref. [7] is based on the QCD-improved parton model implemented in the 
modified minimal-subtraction (MS) renormalization and factorization scheme in its pure 
form with rif — 5 massless quark flavors, which is also known the as the massless scheme 
[9] or zero-mass variable-flavor-number scheme. In this scheme, the masses m c and m b 
of the charm and bottom quarks are neglected, except in the initial conditions of their 
FF's. This is a reasonable approximation for center-of-mass (cm.) energies \fs ^> m c , m b 
in e + e~ annihilation or transverse momenta pr ^> m c , m b in ep and pp scattering, if 
the respective FF's are used as inputs for the calculation of the cross sections for these 
reactions. Hence, we describe the c, b — > X c transitions by nonperturbative FF's, as is 
usually done for the fragmentation of the up, down, and strange quarks into light hadrons. 

The outline of this paper is as follows. In Sec. 2, we briefly recall the theoretical 
framework underlying the extraction of FF's from the e + e~ data, which has already been 
introduced in Refs. [7,10]. In Sec. 3, we present the D°, D + , Df, and A+ FF's we obtained 
by fitting the respective LEP1 data samples from OPAL [8] at leading order (LO) and next- 
to-leading order (NLO) in the massless scheme and discuss their properties. In Sec. 4, we 
present predictions for the inclusive production of these X c hadrons in nonresonant e + e~ 
annihilation at lower cm. energies and compare them with data from other experiments. 
Our conclusions are summarized in Sec. 5. 
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2 Theoretical Framework 



Our procedure to construct LO and NLO sets of D FF's has already been described in 
Refs. [7,10]. As experimental input, we use the LEP1 data from OPAL [8]. 

In e + e~ annihilation at the Z-boson resonance, X c hadrons are produced either directly 
through the hadronization of charm quarks produced by Z — > cc or via the weak decays 
of Xb hadrons from Z — > bb. In order to disentangle these two production modes, the 
authors of Ref. [8] utilized the apparent decay length distributions and energy spectra of 
the X c hadrons. Because of the relatively long X^-hadron lifetimes and the hard b — > Xb 
fragmentation, X c hadrons originating from X b -hadron decays have significantly longer 
apparent decay lengths than those from primary production. In addition, the energy 
spectrum of X c hadrons originating from X^-hadron decays is much softer than that due 
to primary charm production. 

The experimental cross sections [8] were presented as distributions differential in 
x = 2E(X c )/y/s : where E(X C ) is the measured energy of the X c -hadron candidate, and 
normalized to the total number of hadronic Z-boson decays. Besides the total X c yield, 
which receives contributions from Z — > cc and Z — > bb decays as well as from light-quark 
and gluon fragmentation, the OPAL Collaboration separately specified results for X c 
hadrons from tagged Z — > bb events. As already mentioned above, the contribution due 
to charm-quark fragmentation is peaked at large x, whereas the one due to bottom-quark 
fragmentation has its maximum at small x. 

For the fits, we use the x bins in the interval [0.15, 1.0] and integrate the theoretical 
cross sections over the bin widths used in the experimental analysis. For each of the 
four charmed-hadron species considered here, X c = D°, D + , Df, A+, we sum over the two 
charge-conjugate states as was done in Ref. [8]. As a consequence, there is no difference 
between the FF's of a given quark and its antiquark. As in Refs. [7,10], we take the 
starting scales for the X c FF's of the gluon and the u, d, s, and c quarks and antiquarks 
to be fio = 2iTi c , while we take fiQ = Im^ for the FF's of the bottom quark and antiquark. 
The FF's of the gluon and the first three flavors are assumed to be zero at their starting 
scale. At larger scales //, these FF's are generated through the usual Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) [11] evolution at LO or NLO. The FF's of the first three 
quarks and antiquarks coincide with each other at all scales \x. 

We employ two different forms for the parameterizations of the charm- and bottom- 
quark FF's at their respective starting scales. In the case of charm, we use the distribution 
of Peterson et al. [12], 



which is frequently used for the FF's of light hadrons. Equation (1) is particularly suitable 
for FF's that peak at large values of x, as is typically the case for c — > X c transitions. 



D c {x,nl)=N 




(1) 



[{l-xf + exf 



In the case of bottom, we adopt the ansatz 



D b (x,t4) = Nx a (l-x)P, 



(2) 
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Since the b — > X c FF is a convolution of the 6 — > X b fragmentation and the subsequent 
X5 — > X c + X decay, it has its maximum at small x values. Therefore, Eq. (1) is less 
suitable in this case. We apply Eqs. (1) and (2) for the FF's of all four X c -hadron species 
considered here. 

The calculation of the cross section (l/a tot )da/dx for e + e~ — > 7/Z — > X c + X is 
performed as described in Ref. [7], in the pure MS subtraction scheme, i.e., without 
the subtraction terms dQ a (x) specified in Eq. (2) of Ref. [13]. All relevant formulas and 
references may be found in Ref. [10]. As for the asymptotic scale parameter for five 
active quark flavors, we adopt the LO (NLO) value = 108 MeV (227 MeV) from our 

study of inclusive charged-pion and -kaon production [14]. The particular choice of A^ 
is not essential, since other values can easily accommodated by slight shifts of the other 
fit parameters. As in Refs. [7,10], we take the charm- and bottom-quark masses to be 
m c — 1.5 GeV and m b = 5 GeV, respectively. 



3 Determination of the D°, £>+, £>+, and A+ FF's 

The OPAL Collaboration [8] presented x distributions for their full D°, D + , Df, and A+ 
samples and for their Z bb subsamples. We received these data in numerical form via 
private communication [15]. They are displayed in Figs. 4 (for the D° and D + mesons) 
and 5 (for the Df meson and the A+ baryon) of Ref. [8] in the form (l/N h3id )dN/dx : 
where N is the number of A c -hadron candidates reconstructed through appropriate decay 
chains. In order to convert this into the cross sections (l/cr tot )da/dx, we need to divide 
by the branching fractions of the decays that were used in Ref. [8] for the reconstruction 
of the various X c hadrons, namely, 

B(D°^K-7i + ) = (3.84 ±0.13)%, 
B(D + -> R-Ti+Ti-) = (9.1 ±0.6)%, 
fi(D+^07T + ) = (3.5 ±0.4)%, 

B(k+ ^ P K~n + ) = (4.4 ±0.6)%, (3) 

respectively. The experimental errors on these branching fractions are not included in our 
analysis. 

The values of N and e in Eq. (1) and of N, a, and (3 in Eq. (2) which result from our 
LO and NLO fits to the OPAL data are collected in Table 1. From there, we observe that 
the parameters a and /3, which characterize the shape of the bottom FF, take very similar 
values for the various X c hadrons, which are also similar to those for the D* + meson listed 
in Table I of Ref. [7] . On the other hand, the values of the e parameter, which determines 
the shape of the charm FF, significantly differ from particle species to particle species. In 
the D*+ case [7], our LO (NLO) fits to ALEPH [5] and OPAL [6] data, which required 
separate analyses, yielded e = 0.144 (0.185) and 0.0851 (0.116), respectively. We observe 
that, for each of the X c -hadron species considered, the LO results for e are considerably 
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Table 1: Fit parameters of the charm- and bottom-quark FF's for the various X c hadrons 
at LO and NLO. The corresponding starting scales are /io = 2m c = 3 GeV and fiQ = 
2mb = 10 GeV, resp ectively. All other FF's are taken to be zero at /i p = 2m c . 





Order 


Q 


N 


a 




e 


D° 


LO 


c 


0.998 


— 


— 


0.163 






b 


71.8 


1.65 


5.19 


— 




NLO 


c 


1.16 


— 


— 


0.203 






b 


97.5 


1.71 


5.88 





D + 


LO 


c 


0.340 






0.148 






b 


48.5 


2.16 


5.38 






NLO 


c 


0.398 






0.187 






b 


64.9 


2.20 


6.04 




Df 


LO 


c 


0.0704 






0.0578 






b 


40.0 


2.05 


4.93 






NLO 


c 


0.0888 






0.0854 






b 


21.8 


1.64 


4.71 






LO 


c 


0.0118 






0.0115 






b 


44.1 


1.97 


6.33 






NLO 


c 


0.0175 






0.0218 






b 


27.3 


1.66 


6.24 





smaller than the NLO ones. Furthermore, we notice a tendency for the value of e to 
decrease as the mass {m Xc ) of the X c hadron increases. 

In Table 2, we list three values of x 2 P er degree of freedom (xdf) f° r each of the fits 
from Table 1: one for the Z —> bb subsample, one for the total sample (sum of tagged- 
cc, tagged-66, and gluon-splitting events), and an average one evaluated by taking into 
account the Z — > bb subsample and the total sample. The actual Xdf values are rather 
small. This is due to the sizeable errors and the rather limited number of data points, 
especially for the and A+ data. In each case, the Z — > bb subsample is somewhat less 
well described than the total sample. The NLO fits yield smaller Xdf values than the LO 
ones, except for the case. 

The normalized differential cross sections (l/cr tot )da/dx for D°, D + , Df, and A+ 
hadrons (circles), extracted from Ref. [8] as explained above, are compared with our LO 
(upmost dashed lines) and NLO (upmost solid lines) fits in Figs. l(a)-(d), respectively. 
The same is also done for the Z — > bb subsamples (squares). In addition, our LO and 
NLO fit results for the Z — > cc contributions are shown. In each case, the X c hadron 
and its charge-conjugate partner are summed over. From Figs. l(a)-(d), we observe that 
the LO and NLO results are very similar, except for very small values of x. This is 
also true for the distributions at the starting scales, as may be seen by comparing the 
corresponding LO and NLO parameters in Table 1. The branching of the LO and NLO 
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Table 2: x 2 per degree of freedom achieved in the LO and NLO fits to the OPAL [8] data 
on the various D hadrons. In each case, Xdf i s calculated for the Z — > bb sample (6), the 
full sample (All), and the combination of both (Average). 





Order 


b 


All 


Average 


D° 


LO 


1.16 


0.688 


0.924 




NLO 


0.988 


0.669 


0.829 


D + 


LO 


0.787 


0.540 


0.663 




NLO 


0.703 


0.464 


0.584 


Dt 


LO 


0.434 


0.111 


0.273 




NLO 


0.348 


0.108 


0.228 


A c + 


LO 


1.05 


0.106 


0.577 




NLO 


1.05 


0.118 


0.582 



results at small values of x indicates that, in this region, the perturbative treatment ceases 
to be valid. This is related to the phase-space boundary for the production of X c hadrons 
at x m i n = 2m x r J y/s- These values are somewhat larger than the x values where our NLO 
results turn negative. Since our massless- quark approach is not expected to be valid in 
regions of phase space where finite-mx c effects are important, our results should only 
be considered meaningful for x^x CVLt = 0.1, say. We also encountered a similar small-a; 
behavior for the D*+ FF's in Refs. [7,10]. 

As mentioned above, we take the FF's of the partons g, u, u, d, d, s, s to be vanishing at 
their starting scale /iq = 2m c . However, these FF's are generated via the DGLAP evolu- 
tion to the high scale ji = \/s. Thus, apart from the FF's of the heavy quarks c, c, b, b, also 
these radiatively generated FF's contribute to the cross section. All these contributions 
are properly included in the total result for (l/a tot )da/dx shown in Figs. l(a)-(d). At 
LEP1 energies, the contribution from the first three quark flavors is still negligible; it is 
concentrated at small values of x and only amounts to a few percent of the integrated 
cross section. However, the contribution from the gluon FF, which appears at NLO in 
connection with qqg final states, is numerically significant. As in our previous works 
[7,10], motivated by the decomposition of (l/a tot )da/dx in terms of parton-level cross 
sections, we distributed this contribution over the Z cc and Z — > bb channels in the 
ratio e 2 c : ef, where e q is the effective electroweak coupling of the quark q to the Z boson 
and the photon including propagator adjustments. This procedure should approximately 
produce the quantities that are compared with the OPAL data [8]. 

As in Refs. [7,10], we study the branching fractions for the transitions 
c,b^D°,D+, D+, A+ defined by 

£? Q (/i)= f 1 dxD Q (x,pi 2 ), (4) 

where Q = c,b, Dq are the appropriate FF's, and x cut = 0.1. This allows us to test the 
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Table 3: Branching fractions (in %) of c, b — > D°, D + , Df, A+ evaluated according to 
Eq. (4) in LO and NLO at the respective starting scales /i = 2m Q and at the Z-boson 
resonance \i = M z . 



X c Order B c (2m c ) B C {M Z ) B b (2m b ) B b (M z ) 





LO 


72.1 


66.9 


57.8 


52.8 




NLO 


69.5 


63.9 


55.2 


49.8 


D + 


LO 


26.6 


24.7 


19.4 


17.9 




NLO 


25.6 


23.6 


18.6 


17.1 


Df 


LO 


11.5 


10.9 


22.4 


20.6 




NLO 


10.8 


10.1 


21.6 


19.6 


A+ 


LO 


5.88 


5.67 


15.1 


13.7 




NLO 


5.74 


5.48 


14.5 


13.0 



consistency of our fits with information presented in the experimental paper [8] that was 
used for our fits. The contribution from the omitted region < x < x cut is small. Table 3 
contains the values of Bq({i) for all eight transitions c,b — > D°,D + ,Df,Af evaluated 
according to Eq. (4) in LO and NLO at the respective thresholds /x = 2m Q and at the 
Z-boson resonance \x = M z . As expected, the values of Bq(^l) change very little under 
the evolution from \x = Iraq to fi = M z , and they are rather similar for Q = c,b. Leaving 
aside the insignificant contribution due to strange charm baryons, the values of Bq(h) for 
X c = D°, D + , Df, A+ should approximately add up to unity for each heavy flavor Q = c,b 
at any value of /i. Although we did not impose this sum rule as a constraint on our fits, 
it is well satisfied for B C (M Z ) and B b (M z ) at NLO. In fact, from Table 3 one obtains 
103% and 99.5%, respectively. The corresponding LO values, being 108% and 105%, are 
somewhat too large, as may be understood by observing the excess of the LO fits over 
the experimental data at small values of x in Figs. l(a)-(d). The corresponding sums of 
the LO and NLO entries for B c {2m c ) and B b (2m b ) in Table 3 range between 110% and 
116%. In view of the long evolution paths from the charm and bottom thresholds way up 
to the Z-boson resonance, such violations of the sum rule can be considered acceptable. 
The situation is expected to improve once experimental data at lower cm. energies (see 
Sec. 4) are included in our fits. 

It is interesting to compare our LO and NLO values of B C {M Z ) and B b (M z ) for 
the D° ,D + , Df, and A+ hadrons with the respective results determined by the OPAL 
Collaboration through Peterson model fits. These results are presented in Table 9 of 
Ref. [8] in the dressed form 

p Q ^ Xc = RqBq(M z )B Xc , (5) 

where Rq = Tqq /Thad are the production rates of the quarks Q = c, b in e + e _ an- 
nihilation on the Z-boson resonance and Bx c are the decay branching fractions of the 
four X c hadrons X c = D°,D + ,Df,Af into the channels considered in Eq. (4). For 
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Table 4: X c -hadron production rates reported by OPAL [8] compared to results evaluated 
at LO and NLO from Eq. (5) using the branching fractions from Table 3. 



x c 


Pc^X c [%] 

[8] 


fit LO 

nt NT.n 


Pb^x c [%] 
[8] 


fit LO 
nr NT,n 


n(Z -> X c )B Xc 
[8] 


[%] 

fit LO 
nt NT.n 


D° 


0.389 ±0.027 to oil 


0.434 
0.414 


0.454 ±0.023+^26 


0.439 
0.414 


1.784 ±0.066 ±0.086 


1.746 
1.656 


D+ 


0.358 ± 0.046 tooti 


0.380 
0.363 


0.379 ± 0.031 +°'o25 


0.353 
0.337 


1.548 ±0.082+°'°^ 


1.466 
1.400 


D+ 


0.056 ±0.015 ±0.007 


0.0644 
0.0597 


0.166 ±0.018 ±0.016 


0.156 
0.148 


0.460 ± 0.036 ± 0.040 


0.441 
0.415 


A+ 


0.041 ±0.019 ±0.007 


0.0421 
0.0407 


0.122 ±0.023 ±0.010 


0.130 
0.124 


0.345 ±0.052 ±0.029 


0.344 
0.329 



the reader's convenience, these results are copied to Table 4, where they are compared 
with our results for pq^ Xc , which are obtained from the appropriate entries in Table 3 
through multiplication with the branching fractions from Eq. (4) and the production rates 
R c = 0.1689 ±0.0047 and R b = 0.21643 ±0.00072 determined by the Particle Data Group 
[16] in the framework of the Standard Model. For simplicity, the values deduced from 
Table 4 do not include the errors on Rq and Bx c and those on Bq(M z ) resulting from 
our fits. 

In Table 9 of Ref. [8], the OPAL Collaboration also presented the total rates n(Z — > 
X c )B Xc , which include the estimated contributions from gluon splitting g — > QQ; for 
further details, see Ref. [8]. In Table 4, these results are quoted and compared with 
the corresponding quantities 2(p c ^x c ± Pb~*x c ) resulting from our LO and NLO analyses. 
Notice that the experimental results are corrected to include the unmeasured contributions 
from x < 0.15, whereas our evaluations of Eq. (4) exclude the contributions from x < x cut . 
This explains why the experimental results somewhat overshoot ours. The agreement is 
worse at NLO, which may be understood by observing that our evaluations of 2{p c ^x c + 
Pb^x c ) do not include the contributions from gluon fragmentation, which enters the stage 
at NLO. Keeping these caveats in mind, we find reasonable overall agreement between 
the OPAL results for n(Z — > X c )B Xc and our results for 2(p c _^ Xc + Pb-+x c )- 

Our LO and NLO values of B C {M Z ) and B b (M z ) for the D°, D + , D+, and A+ hadrons 
in Table 3 can also be compared with experimental results published more recently by the 
ALEPH [5] and DELPHI [17] Collaborations. In Ref. [5], B C (M Z ) are called f(c -> X c ) 
and may be found in Sees. 7.1 and 7.3. In Ref. [17], Bq(M z ) are called Pq^x c and may 
be extracted for Q = c from Table 13 (in connection with sum rule of Eq. (12) and taking 
into account the discussion of the contribution from the strange charm baryons in Sec. 8.2) 
and for Q = b from Table 15. For simplicity, we add the three types of errors quoted in 
Refs. [5,17] (from statistics, systematics, and decay branching fractions) in quadrature. 
In 1999, Gladilin [18] derived world-average values of B C (M Z ) for the D° ,D + , D+, and 
A+ hadrons related to e + e _ annihilation, which are also listed in Table 5. 

The branching fractions of the c — > D°,D + ,Df,A+ transitions were also measured 
in ep collisions at HERA, in photoproduction by the ZEUS Collaboration [1,2] and in 
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Table 5: Branching fractions (in %) of c, b -> Z>°, L>+, D+, A+ reported by ALEPH [5], 
D ELPHI [17], Gladilin [18], ZEUS [2], and HI [3]. 

5 C (M Z ) B b {M z ) 

[5] [17] [18] [2] [3] [17] 

~D° 55.9 ±2.2 54.80 ±4.78 54.9 ± 2.6 55/r^ 65.8™ 60.05 ±4.39 

L>+ 23.79 ±2.42 22.70 ± 1.82 23.2 ± 1.8 24.9+j'^ 20.2 "I 23.01 ±2.13 

£>+ 11.6 ±3.6 12.51 ±2.97 10.1 ±2.7 10.7 ±1.0 15.6+^ 16.65 ±4.50 

A+ 7.9 ±2.2 8.76 ±3.30 7.6 ±2.1 7.6+^g - ' 8.90 ± 3.00 



deep- inelastic scattering by the HI Collaboration [3]. These results are also included in 
Table 5 for comparison. Strictly speaking, they do not correspond to B c (Mz), but rather 
to B c (n), where \x is set by the average value of px (in the case of photoproduction) or 
Q (in the case of deep-inelastic scattering). However, from Table 3 we know that the \i 
dependence of B c {p) is relatively mild. 

We observe that the experimental results collected in Table 5, which are mostly in- 
dependent from each other, are mutually consistent within errors. Comparing them with 
the corresponding entries in the forth and sixth columns of Table 3, we find resonable 
overall agreement. 

Another quantity of interest, which can directly be compared with experiment, is the 
mean momentum fraction, 

1 f 1 

(x) Q {n) = — — - / dxxD Q {x,n). (6) 

In Table 6, we present the values of (x)q(/j,) for Q = c,b evaluated at \x = 2itiq, Mz with 
the LO and NLO FF's of the D°, D + , D+, and A+ hadrons. At fixed value of fi, the 
differences between the LO and NLO sets are insignificant. The DGLAP evolution from 
\x = 2rriQ to /i — M z leads to a significant reduction of (x)q(ii), especially in the case 
of Q = c. The values of (x)b(ii) are appreciably smaller than the values of (x) c (/i), as is 
expected because the bottom-quark fragmentation into X c hadrons is much softer than 
the charm-quark one. 

Our values of (x) c (M z ) for the D° and D + mesons should be compared with the re- 
spective results obtained by the OPAL Collaboration [8] in the framework of the Peterson 
model [12], which read 

(D+) (7) 
for the D° and D + mesons, respectively. The differences to the values obtained for three 



(x) c (Mz) 
(x) c (M z ) 



0.487 ±0.009 
0.483 ±0.015 



±0.011 

-0.009 
±0.007 

-0.011 
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Table 6: Average momentum fractions of c, b — > D°,D + ,Df,Af evaluated according to 
Eq. (7) in LO and NLO at the respective starting scales \x = 2m Q and at the Z-boson 
resonance \i = Mz- 



X c Order (x) c (2m c ) (x) c (M z ) (x) b (2m b ) (x) b (M z ) 





LO 


0.588 


0.452 


0.316 


0.284 




NLO 


0.568 


0.431 


0.300 


0.270 


D + 


LO 


0.596 


0.458 


0.341 


0.303 




NLO 


0.575 


0.436 


0.323 


0.287 


Df 


LO 


0.676 


0.512 


0.349 


0.310 




NLO 


0.644 


0.482 


0.332 


0.296 


A+ 


LO 


0.791 


0.590 


0.302 


0.273 




NLO 


0.750 


0.553 


0.288 


0.261 



other fragmentation models are included in the systematical errors. Comparing Eq. (7) 
with the corresponding entries in Table 6, we observe that the latter are slightly smaller. 

4 Comparison with e + e~ data at lower energies 

The fractional energy spectra of inclusive D°, D + , Df, and A+ production was also 
measured in nonresonant e + e~ annihilation at lower energies. Specifically, the CLEO 
Collaboration took D°, D+ (Table XII in Ref. [19]), Df (Table IV in Ref. [20]), and A+ 
(Table V in Ref. [21]) data at LEPP CESR with = 10.55 GeV; the HRS Collaboration 
took D°, D+ (Table 1 in Ref. [22]), and Df (Table I in Ref. [23]) data at SLAC PEP 
with = 29 GeV; and the TASSO Collaboration took Df (Fig. 3 in Ref. [24]) data at 
DESY PETRA with y/s = 34.7 GeV. It is instructing to confront these data with LO and 
NLO predictions based on our new FF's, so as to test the scaling violations predicted by 
the DGLAP evolution equations. An especially interesting situation arises for the CLEO 
data [19-21], from which all X c hadrons coming from Xft-hadron decays are excluded by 
appropriate acceptance cuts, so that only rif = 4 quark flavors are active and a direct test 
of the charm-quark FF's is feasible. 

The D + and A+ data explicitly refer to the decay channels Df — > 07r + — > K + K~n + 
and A^ — > pK~7r + , respectively, and we have to divide them by the corresponding branch- 
ing fractions. For this, we use the up-to-date values B{Df — > <j)ir + )B((j) — > K + K~) = 
(3.6±0.9)%and5(A+ -> pK-n + ) = (5.0±1.3)% [16], except for the Df data of Ref. [24]. 
In the latter case, for consistency, we adopt the value B(Df — > 07r + )_B(0 — > K + K~) = 
0.13 ± 0.03 ± 0.04 from Ref. [24] itself. 

The differential cross sections da/dx for the D°, D + , Df, and A+ hadrons measured 
by CLEO [19-21] (circles), HRS [22,23] (squares), and TASSO [24] (diamonds) are con- 
fronted with our LO (dashed lines) and NLO (solid lines) predictions in Figs. 2(a)-(d), 
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respectively. 

Let us first concentrate on our NLO predictions. As for the D°, D + , and Df mesons, 
we observe that our NLO predictions generally lead to a satisfactory description of the 
experimental data, both in normalization and shape. In particular, the maxima of the 
measured x distributions are approximately reproduced. However, in the case of the 
A+ baryon, the predicted x distribution appears to be too hard, its peak being set off by 
approximately +0.2 relative to the one shown by the experimental data. In particular, the 
data points at 0.55 and 0.65 are poorly described by the NLO prediction. Although the A+ 
baryon is 22% heavier than the D° and D + mesons, and 16% heavier than the meson, 
mass effects are unlikely to be responsible for this disagreement, since y/s = 10.55 GeV 
is sufficiently far above the charm threshold. 

Let us now include the LO predictions in our considerations. The CLEO data [19-21], 
which are most precise, clearly favor the NLO predictions, while the LO predictions are 
too large at small values of x and too small in the peak region. Unfortunately, the HRS 
[22,23] and TASSO [24] data do not reach the small-x regime, where the LO and NLO 
predictions depart from each other, and their errors are too large in order to support this 
observation. 

Actually, the CLEO data [19-21] are considerably more precise than the OPAL data 
[8], which we fitted to, and it would be desirable to also include them in ours fits. However, 
we refrain from doing so for the time being because their high precision would make it nec- 
essary to properly treat finite-mQ effects, which are neglected altogether in the theoretical 
formalism employed here. The general-mass variable-flavor-number (GM-VFN) scheme 
[25], which has recently been extended to inclusive XQ-hadron production in 77 [26], ep 
[27], and pp [28] collisions, provides a rigorous theoretical framework that retains the full 
finite-mQ effects while preserving the indispensible virtues of the factorization theorem 
[29], namely the universality and the DGLAP [11] scaling violations of the FF's entailing 
the resummation of dominant logarithmic corrections. A global analysis of experimental 
data on inclusive X c -hadron production in the GM-VFN scheme is left for future work. 

5 Conclusions 

The OPAL Collaboration presented measurements of the fractional energy spectra of 
inclusive D°, D + , D+, and A+ production in Z-boson decays based on their entire LEP1 
data sample [8] . Apart from the full cross sections, they also determined the contributions 
arising from Z — > bb decays. This enabled us to determine LO and NLO sets of FF's for 
these X c hadrons. 

As in our previous analysis of D* + FF's [7], we worked in the QCD-improved parton 
model implemented in the pure MS renormalization and factorization scheme with nf = 5 
massless quark flavors (zero-mass variable-flavor-number scheme). This scheme is partic- 
ularly appropriate if the characteristic energy scale of the considered production process, 
i.e., the cm. energy y/s in the case of e + e~ annihilation and the transverse momentum 
Pt of the X c hadron in other scattering processes, is large compared to the bottom-quark 
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mass m b . Owing to the factorization theorem, the FF's defined in this scheme satisfy 
two desirable properties: (i) their scaling violations are ruled by the timelike DGLAP 
equations; and (ii) they are universal. Thus, this formalism is predictive and suitable for 
global data analyses. 

We verified that the values of the branching and average momentum fractions of the 
various c, b — > X c transitions evaluated at LO and NLO using our FF's are in reason- 
able agreement with the corresponding results from OPAL [8] and other experiments 
[2,3,5,17,18]. 

We tested the scaling violations of our FF's by comparing the fractional energy spectra 
of inclusive D°, D + , D+, and A+ production measured in nonresonant e + e~ annihilation at 
yfe = 10.55 GeV [19-21], 29 GeV [22,23], and 34.7 [24] with our LO and NLO predictions 
to find reasonable agreement. Since events of X c -hadron production from X^-hadron 
decay were excluded from the data samples at ^Js = 10.55 GeV, we obtained a clean test 
of our charm-quark FF's. 

It is important to bear in mind that the fit results for the input parameters in Eqs. (1) 
and (2), including the value of Peterson's e parameter, are highly scheme dependent at 
NLO, and must not be naively compared without careful reference to the theoretical 
framework which they refer to. 
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Figure 1: The normalized differential cross sections (l/a tot )da/dx of inclusive (a) D°/D , 
(b) D ± 1 (c) Df, and (d) production in e + e~ annihilation on the Z-boson resonance 
evaluated at LO (dashed lines) and NLO (solid lines) with our respective FF sets are 
compared with the OPAL data [8] renormalized as explained in the text (circles). The 
same is also done for the Z — > bb subsamples (squares). In addition, our LO and NLO 
fit results for the Z — > cc contributions are shown. In each case, the X c hadron and its 
charge-conjugate counterpart are summed q^er. 
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Fig. 1 (continued). 
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Figure 2: The differential cross sections da/dx (in nb) of inclusive (a) D°/D , (b) D ± , 
(c) Df, and (d) Af production in e + e~ annihilation at yHi = 10.55, 29, and 34.7 GeV 
evaluated at LO (dashed lines) and NLO (solid lines) with our respective FF sets are 
compared with data from CLEO at CESR (squares), HRS at PEP (circles), and TASSO 
at PETRA (diamonds), respectively. In each case, the X c hadron and its charge-conjugate 
counterpart are summed over. 
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